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Abstract 

We report the direct observation of a peculiar lava channel that was formed near the base of 
a parasitic cone during the 2001 eruption on Mount Etna. Erosive processes by flowing lava are 
commonly attributed to thermal erosion. However, field evidence strongly suggests that models of 
thermal erosion cannot explain the formation of this channel. Here, we put forward the idea that 
the essential erosion mechanism was abrasive wear. By applying a simple model from tribology 
we demonstrate that the available data agree favorably with our hypothesis. Consequently, we 
propose that erosional processes resembling the wear phenomena in glacial erosion are possible in 
a volcanic environment. 

PACS numbers: 46.55.+d, 81.40.Pq, 91.40.Hw, 92.40. Gc 
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Introduction - Formation of lava channels is usually discussed in terms of thermal or 
thermomechanical erosion DDQDDQD . The basic idea is that heat is transferred from 
the hot flowing lava to the cold substrate, thereby melting it and washing it out. Thermal 
erosion by lava has been studied quantitatively both for turbulent and for laminar 

flows . The physics underlying the quantitative description is that of heat transport in 

liquids and solids combined with fluid dynamics. The idea of thermal erosion in the formation 
of lava channels is supportec 
in archaean komatiite lavas 
Kilauea volcano (Hawaii) 

In the picture of thermal erosion, the relevance of mechanical processes appears to be 
secondary being mainly limited to the removal of the melted material. However, it has been 
noted that mechanical effects may enhance the erosion rates 0, Q] . Very recently, Williams 
et al. have reported geochemical evidence for mechanical erosion in the Cave Basalt of 
Mount St. Helens where a basaltic flow eroded an unconsolidated pyroclastic substrate jl^j ]. 
To our knowledge, this is the only example in the literature in which a dominant role of 
mechanical erosion of lava has been inferred (aspects of models for mechanical erosion have 

). Therefore, the direct observation of the formation 
of a lava channel whose origin cannot be explained by thermal erosion is of considerable 
interest. Such an observation may help to develop new and more complete models for the 
explanation of the field data elsewhere. Moreover, it generates interesting cross-links to 
other branches of geology where a variety of erosion models are discussed. 

The starting point of this work was the formation of an uncommon channel during the 
2001 eruption on Mount Etna which provides field evidence for purely mechanical erosion 
of lava flow. First, we briefly describe the morphological details of the channel and the 
history of the eruption (see also Ref. An estimate of the heat transfer required for the 

formation of the channel basically rules out thermal erosion as the essential process. Then, 
we will discuss the possibility of erosion due to abrasive wear and demonstrate that this 
hypothesis may account also quantitatively for the field observations. We will conclude the 
discussion of our model by deriving some theoretical predictions and compare them to those 
for thermal erosion. 

Description of the Laghetto channel - The channel, informally named Laghetto, is located 
on the upper southern flank of Mount Etna volcano at an elevation of 2560 m a.s.l. (see Fig. 
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1). It is roughly linear, of rectangular-shaped cross section and is incised on a lava slope 
gently inclined (5°-10°). Its axis is oriented E-W. It is 220 m long and has a width between 
8 and 16 m. Maximum depth is 6 m, the floor of the channel is covered by lava. The banks 
consist of lava flows emitted between July 25th and 30th. They are formed of up to 11 thin 
lava layers, each about 0.2 to 0.3 m thick (cf. Fig. 2). The eroded volume is about 220 m x 
12 m x 4 m k 10,000 m 3 . Note that this volume is underestimated since part of the channel 
is filled with the eroding lava. 

Description of the eruptive event - The eruption started on July 19th forming a crater in 
the area of Piano del Lago. In the first week, the eruption (closely observed by the authors) 
gave rise only to phreatomagmatic explosions and lava fountains. From July 25th for the 
5 following days lava was erupted from the vent. The relatively low lava output (probably 
< 5 m 3 /s) was not constant but intermittent thus originating a number of thin and short 
(average length < 500 m) lava flows. 

In the evening of July 30th we witnessed an abrupt change in the lava output. A huge 
lava flow (10 to 20 m high) was emitted. The lava flowed (in laminar motion) for about 12 
hours reaching the area of Rifugio Sapienza, 2000 m down valley. Near the vent the lava 
flowed on a bedrock formed by the still hot and soft thin lava layers of the preceding days. 
In these hours the bedrock was eroded and the following morning the Laghetto channel was 
there. After this last flow, lava emission at the Laghetto vent ceased. Therefore, the channel 
was not buried under subsequent lavas. It is this fortunate circumstance that allowed us to 
study the phenomenon and measure the channel before it was partially covered by sand and 
ash fall of the flank eruption of October 2002. 

Thermal erosion hypothesis - Let us assume that the channel was thermally eroded. This 
implies that the lava flow had to supply (at least) the energy A£ to melt the material 
originally contained in the channel, i.e., the volume W v ~ 10 4 m 3 

AS = pW v ■ ( L + c(T m — T subst )) . (1) 

Here c = 1.4 ■ 10 3 J/ (kg K) is the specific heat of the lava substrate (calculated for the 
specific lava composition using MELTS ^f|), T m « 1050°C is a lower limit for the onset of 
melting, and T su b s t ~ 900° C denotes the likely substrate temperature. With a lava density of 
p = 2800 kg/m 3 and a latent heat of L = 0.5 • 5.0 • 10 5 J/kg for an assumed 50% cristallinity 
(using typical values for basaltic lavas 0, Q, Q), about 1.3 x 10 13 Joule were required. 
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The eroding lava flow had a total length s ~ 2000 m (the distance between vent and 
the Rifugio Sapienza area) and an average height ha ~ 10 m. Typically, Mount Etna lavas 
are rich in phenocrysts (i.e., they are erupted partially solidified) and their temperature 
(T, ~ 1050 - 1100°C) is close to the solidus point. 

Assuming an excessive temperature difference 7] — T m ~ 50 K, the lava flow (approx. 
240,000 m 3 ) could provide at most 4.7 x 10 13 Joule of heat energy (during 12 hours). There- 
fore, a significant part of the available heat energy had to be transferred extremely rapidly 
from the lava flow to the substrate to thermally erode the channel. The corresponding 
heat current density (transferred energy/ (time x channel surface)) amounts to ~ 10 5 W/m 2 ). 
Even with a temperature gradient on the order of ~ 100 K/m this corresponds to a heat 
conductivity of about ~ 10 3 W/(m-K). 

This is rather unrealistic as can be seen by the following considerations. First, the huge 
heat transfer rate would have required a heat conductivity exceeding by far that of noble 
metals {s-SLj, f° r silver one has n\ g 4 x 10 2 W/(m-K) [12( while basaltic lavas have k ~ 1 
W/(m-K) jy, 2|). Second, the noticeable heat loss would have increased viscosity leading to 
rapid deceleration and eventual standstill of the flow. This is in contradiction with the field 
observations. Finally, we have neglected heat loss, e.g., radiative heat loss and the formation 
of chill layers Thus, we conclude that thermal erosion cannot have played a dominant 
role in the channel formation. 

Mechanical erosion hypothesis - Friction and wear are ubiquitous phenomena in every- 
day life. Their investigation in physics and engineering sciences, however, is a relatively 
young field Q, Q. Mechanical erosion is one of the fundamental processes in geology, 



21] erosion. While 



quantitatively studied, e.g., in the context of bedrock [2Cj and glacial 
such studies focus on the dynamics, based on field data for many examples, the aim of the 
present discussion is merely to argue that the formation of the Laghetto channel does comply 
with the fundamental concepts of abrasive wear. 

The idea of purely mechanical erosion is imposed by the fact that strata in the channel 
banks are clearly distinguished (cf. Fig 2), indicating that the material was removed by 
"cutting" rather than by melting. The question is which mechanisms have acted and how 
they can be described quantitatively. 

We propose the following intuitive picture. The erosive lava was a thick (up to 20 m) 
autoclastic flow, i.e., a dense mixture of crystals and clasts, metric in size, suspended in a 
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melt matrix. The mesh of clasts got dragged along the bed such that the irregular clast 



surfaces in contact with the bed cou 
abrasive and erosive wear processes 



d plough out grooves. These are the characteristics of 



221. 



It needs to be emphasized that abrasion and wear are rather complex phenomena and 
therefore hard to model at a 'microscopic' level (that is, by identifying a dominant mech- 
anism and deducing, e.g., wear rates from dynamical equations and material properties, in 
particular if the process parameters cannot be defined as well as in systematic lab exper- 
iments). For this reason, phenomenological models that are well-established by empirical 

studies and describe whole classes of wear processes are of considerable importance. One 

I I 

of the simplest quantitative models to describe wear phenomena is Archard's law [18]. It 
establishes a phenomenological relation between the wear volume W v , the normal load F, 
and the sliding distance s: 

Here, H denotes the hardness of the worn material proportional to the yield strength r 
(H ~ 3t). Further, k is the wear coefficient, i.e., a dimensionless constant which depends on 



the wear mechanism. Abrasive wear by sliding glaciers is described by similar models [21 1. 
From the observations we can estimate W v , s, F, and H. For a realistic model, k should be 
0.1 or smaller. The wear volume is W v ~ 10, 000 m 3 . The sliding distance is the total length 
of the flow (s ~ 2000 m) and the load F is the weight of the abrading lava along the length 
of the channel. With an average height of the flow of 10 m we find F w 26, 000 m 3 x 2800 
kg/m 3 x 9.8 m/s 2 » 7.3 x 10 8 N. 

It is more difficult to estimate the hardness H of the substrate since its yield strength is 
an exponential function of temperature. Note that these lava layers had been erupted only 
days or hours before. Their outer skin was black, but they were sti 
For cold lavas yield strength values are ~ 10 7 Pa. According to Ref. 
of 1000 °C yield strength is larger than 10 4 Pa. For a likely temperature of 900 °C we 
interpolate values between 10 5 to 10 6 Pa. 

By inserting these numbers into Eq. (J2J) we obtain k « 10~ 3 . . . 10~ 2 . This wear coefficient 
lies well within the typical range for abrasive and erosive wear processes (k ~ 10 -5 . . . 10 _1 ) 



1 incandescent inside. 
23| for a temperature 



221 ]. confirming the intuitive picture for the erosive action of the lava 



measured in tribology 
flow. 

Note that our estimate of k and its relation to the wear mechanism is rather robust 



5 



against variations due to uncertainties in the observations while thermal erosion remains 
excluded by several orders of magnitude. 

In view of the strong support that the experimental observations provide for the idea of 
mechanical erosion, we may invert the argument: we may ask what theoretical predictions 
can be obtained from assuming this kind of erosive mechanism, and compare them with 
the results for thermal erosion. To this end, we rewrite Archard's law Eq. (J2J) in terms 
of characteristic quantities of the erosive process. We have already mentioned the sliding 
distance s, the height of the flow Jiq, the density of the lava p, the hardness H of the substrate, 
and the constant k which characterizes the wear mechanism. Further we introduce the width 
Wa of the flow (the latter is equal to the channel width w c h), the depth d c ^ and the length 
Z c h of the channel. With W v = Z c h^ch^ch and F = pgl^w^h^cosa (where a is the angle of 
the slope; for simplicity we will approximate cos a ~ 1 for moderate slopes) we obtain the 
relation 

d c h = frfl s ■ (3) 

This equation can be further transformed by assuming that the flow rate per unit width Q 
was constant for the duration ig of the flow. The result is 

Knech = Q , (4) 

where V mec h denotes the velocity of mechanical erosion V mec h = .f-rfch- This relation can be 
compared to Kerr's key result for the melting velocity Vthermai due to laminar flows |2j 



U \ 1/3 

Vthermai OC ( J , (5) 

^fl, thermal / 



here U denotes the surface velocity of the flow, and we have identified = ha for thermal 
erosion. While Eq. © implies, e.g., [3] 

^thermal OC Q 1 ^ , /^thermal OC Q 1 ^ 3 (6) 

our simple theory of mechanical erosion predicts 

,mech (7) 

Moreover, the substrate yield strength depends exponentially on the substrate temperature: 
H = aexp (— &T su bst) (with material constants a, b > 0) and thus determines the temperature 
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dependence of the ratio V mcc h/Q in Eq. while for thermal erosion the ratio Vthermai/Q 1//3 
depends on T su b s t via a rational function (the Stefan number, cf. Ref. pj). 

Thus, our hypothesis leads to concrete predictions which are substantially different from 
those for thermal erosion, even without precise knowledge regarding the 'microscopic details' 
of the erosion mechanism. An important consequence of Eqs. (JBJ) and (J7J) is that, with 
increasing flow rate Q, the velocity of mechanical erosion Vm e ch grows much faster than that 
of thermal erosion. That is, for typical parameters of basaltic lava flows ?[ 23 1 one may 
obtain velocities of mechanical erosion that are one or two orders of magnitude larger than 
those for thermal erosion due to laminar flows with the same flow rate. This fact may help to 
explain the observation of lava-channel formation even in the absence of very high eruption 
temperatures and long-lived, low- viscosity flows (cf. Ref. jl^). 

Conclusions - The analysis presented here focused mainly on explaining the physical 
processes at the origin of the Laghetto channel. It may serve as a basic approach to me- 
chanical erosion by lava flows. Further field studies and a better understanding of the erosive 
processes at the 'microscopic level' will help to refine the model. It is important to note 
that the observation of the channel formation described here represents a singular event. 
Nevertheless, similar features indicating mechanical erosion can be found in historical and 
prehistorical lava fields on Mount Etna and elsewhere Finally we mention that new 

aspects in channel formation processes on Earth may provide useful hints also for the as- 
sessment of erosion mechanisms in planetary environments such as lunar sinuous rilles and 
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Fig. 1: Aerial view of the Laghetto eruptive vent. The channel can be seen to the left 
of the pyroclastic cone. As indicated the erosional feature (head of the channel) starts 
immediately at the emission vent. Its length (220 m) gives the scale of the picture. In the 
upper part, there is the lava field produced by the eruption of 2001. 



Fig. 2: The right bank of the channel. In the foreground the bottom of the channel can be 
seen. It is covered by the eroding lava flow. The incision of the channel renders visible the 
layered structure of the substratum. It is important to note that these layers erupted few 
days or hours before were still rather soft and plastic. This is displayed by the uppermost 
layer folded inside the channel just after its incision. 



Fig. 3: Cartoon showing the mechanical erosion of the channel (not in scale), a) Near the 
vent (less than 500 m), the preflow terrain consisted of thin layers of hot lava, b) While 
moving down the slope the autoclastic lava flow carved out the soft parts of the substrate, 
c) Situation after the flow had stopped. The channel floor is partially covered by lava of the 
eroding flow. 
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